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Constrained Mechanics and Variational Problems 
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W'sKO 

A BRIEF SURVEY OF CONSTRAINED MECHANICS AND n > 

VARIATIONAL PROBLEMS IN TERMS OF DIFFERENTIAL FORMS / ' ' 


Robert Her mann 

Ever since my graduate student days, I have been impressed and 
influenced by the elegance and systematization of Mechanics and 
Variational Calculus contained in Ehe Cartan's book "Lecons sur les 
Invariants Integraux”. In the period 1959-69, 1 expended considerable 
effort in the development of Cartan's point of view in many books and 
articles. In this paper (which will appear as a Chapter in ”In terdisci phnay 
Mathematics", v. 30), I will give a quick development of some of the 
material in my books "Differential Geometry and the Calculus of 
Variations" and "Geometry, Physics and Systems”. 

Another purpose in developing this geometric form of the Equations 
of)Mechanics in this Volume is that it fits in with my strategy of 
investigating mechancis with 'singular' features, such as Delta Functions, 
Discontinuities, Shocks, etc. As I will show In Volume 30 the C-O-R 
constructions of Generalized Functions enable one to define 'differential 
forms with generalized coefficients', thus preparing the ground for the 
material in this Chapter serving as foundation for Mechanics with Singular 
Data, the Theory of Splines on nonlinear manifolds, etc. Further, when 
combined with the Computational Methods under development at the AI 
Lab of MIT by Gerry Sussman and co-workers this material will be useful 
in the development of Air Traffic Control methodology. 

Another gaol of my work is to develop a general structure for ODE 
systems, to be used in both 'smooth' and 'gene raliz ed' (in the sense of 
Colombeau, Oberguggenberger and Rosinger) Mechanics, Control and 
Numerical Analysis. Since Martin, Crouch have shown that,-in-the linear 
case, Splines may be constructed from linear control system so attention 
will, in the future, focus on the Splines associated with Generalized Inpouts 
to Nonlinear Con trol Systems. Work ofSastry and Montgomery indicates 
that imporant examples of such systems will be the Left-Invariant Control 
Systems on Lie Groups, which have been much studied in recent years by 
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Constrained Mechanics and Variational Probl 



researchers interested in Integrahle Systems, Robotics, and Aircraft 
Guidance, 



1. Introduction. 

There has been considerable interest recently in constrained 
mechanics and variational problems. This is in part due to Applied 
interests (such as 'non-holonomic mechanics in robotics') and in other part 
due to the fact that several schools of 'pure' mathematics have found that 
this classical subject is of importance for what they are trying to do. I 
have made various attempts [2, 3, 4, 5, 6, 8, 11, 15, 20, 26, 27] at 
developing these subjects since my Lincoln lab days of the late 1950's. In 
this Chapter, I will sketch a Unified point of view, using Cartan's approach 
with differential forms. This has the advantage from the C-O-R viewpoint 
being developed in this Volume that the extension from 'smooth' to 
'generalized' data is very systematic and algebraic. (I will only deal with 
the 'smooth' point of view in this Chapter; I will develop the 'generalized 
function' material at a later point.) The material presented briefly here 
about Variational Calculus and Constrained Mechanics can be found in 
more detail in my boooks, "Differential Geometry and the Calculus of 
Variations" "Lie Algebras and Quantum Mechanics", and "Geometry, Physics 
and Systems". 

Here is the basic set-up. Suppose given the following data: 


A smooth paracompact manifold X (1.1) 

T(X) is its tangent vector bundle (1.2) 

A set {0, co 1 , ..., co m ] of smooth 1-forms on X. (1.3) 

6 is called the action form, {to 1 , ..., co®} are the constraint forms. 


Let us suppose given a curve in X: 
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37 Constrained Mechanics and Variational Problems 

x= {t--> x(t) eX: astsb}}: [a,b] -->X (1.4) 

dx/dt = v = {t~> dx/dt(t) eT(X): astsbj|: [a, b] --> T(X) (1.5) 
is its tangent vector or velocity curve. 

(In this Chapter, I suppose all such curves are also smooth.) 

De fini tion. The following real number associated to the curve 1.3 is called 
the action: 


°(x) =/[a,b]9([dx/dt](t))dt (1.6) 

The following field of 1-covectors along the curve 1.4 is called the force: 

ft > [dx/dt](t))J de} (1.7) 

1.6 and 1.7 are the basic data for both 'mechanics' and 'variational 
calculus'. 

Now, let us deal with 'constraints': 

Definition. The curve 1.4 satisfies the constraints associated with 
the 1-forms {to 1 , ..., a> m } iff. it satisfies the following set of Pfaffian 
differential equations: 

0 = co 1 (dx/dt) = ... = u) m ( dx/dt) (1.8) 

I will show how the basic Equations of Mechanics can be described 
very compactly and elegantly in terms of this data. 

2. The First Variation formula and the Cauchy characteristic 
curves of de. 
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3 8 Constrained Mechanics and Variational Problems 

Keep the data of Section 1. Let us suppose that only the 'action 1 form 
9 is give, without any constraints. Let x be a smooth curve in X, given as in 
1.4. Suppose that 's', 0 Is s s 1, is a deformation parameter and that 
{s — > Xs.* [a, b] — > X} is a smooth one-parameter family of curves in X, 
reducing to the given curve x at T s=0’. For t e [a, b], set: 

v(t) = tangent vector to the curve fs --> Xs(t)} at 't= O’ 

The field {t --> v(t) sX X ( t )} of tangent vectors is called an infini tesimal 
deformation of the curve x. Then: 

d(a(Xs))/ds|s»o (2.1) 

is called the First Variation of the action function function 1.6 
along the curve x pointing in the direction of the vector field v. 
Using the formula 1.6 for the Action, the Cartan Family Identity: 

’V(e) = VJ de + d(VJ 9)’: between a differential form and a vector field, and 
Integration-By-Parts, we have the First Variation Formula: 

d(a(xs))/ds| s =o = /[a, b] -[dx/dtj de](v(t))dt + 0(v)(b) - 0(v)(a) (2.2) 


Remark. This formula is a variant of a General Principle: 

The Variational Derivative of the Action is the Force (2.3) 
It also suggests the following: 

Definition. A curve (x: t --> x(t)| is called a Cauchy characteristic 
curve for the 2 -form d0 iff: 

[dx/dt] J de = 0. " (2.4) 

If x satisfies 2.4, then the First Variation 2.2 vanishes for any infinitesimal 
deformation v which satisifies the following conditions: 
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39 Constrained Mechanics and Variational Problems 

e(v)(b) = e(v)(a) = 0 (2.5) 

Conditions 2.5 are called Transversality Conditions. 

At this point, 'symplectic structures and foliations', 'Hamilton's and 
Lagrange's Equations' (for special choices of 0), etc. enter in a very natural 

way. See [2, 4, 6, 8, 11, 20, 27, 29]. 

3. The differential equations of constrained extrema and the 
augmented action form. 

Let us now suppose that {6, co 1 , ..., to m } is given, as in 1.3. Introduce 
Lagrange Multiplier Variables: 

{Xi, •••» (3.1) 

Consider them as Cartesian coordinates of a copy of R m . On X x R m , 
introduce the following augmented 1-form: 

0aug ~ 0 + XjojI + ... + 

Then, 

d0aug == d0 + dX^Aco^ + ... + dX m A 0 »m +■ 

Xidco 1 + ... + Xmdu/^ 

Definition. A curve {t --> x(t)} in X is an extremal of the constrained 
variational problem associated with the differential form data 1.3 if and 
only if there is a curve in X x R m of the form 

{t -> (x(t), Xi(t), ..., X m (t)} which is a Cauchy characteristic curve of de aug . 
In other words, the 'extremas' are the images under the Cartesian 
projection map: {X x R m -> X] of the Cauchy characteristic curves of de aug . 

Theorem 3.1. A curve {t --> (x(t), Xi(t), ..., Xm(t)} is a Cauchy 
characteristic curve for the 2-form 'd0 au g' if and only if the following 
conditions are satisfied: 


(3.2) 

(3.3) 


0 = to^dx/dt) = ... = (^(dx/dt) 


(3.4) 
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dx/dtj de = - Xi(t)[dx/dt] J do>i . ... . X m (t)[dx/dt] J du>m (3.5) 

- [dxi(t)/dt] col . . [dx m (t)/dt] com 

Proof. Let v be a tangent vector to the manifold X x R m . Then: 

vj dBaug = vj (d0 + dXiAcol + ... + dX m Aco m + Xjdcol + ... + Xjndco 111 ) 

= vj de + v(Xi)o)l + ... + v(Xi)col - col(v)dXi - ... - co m (v)dX m (3.6) 

+ Xi[vJ do 1 ] + ... + X m [vj dco 01 ] 

3.6 involves one-forms on X x R m . Notice that the only terms on the right 
hand side of 3.6 which involves fdXi, ... dx ra } are the terms 
u)!(v)dXi - ... - co m (v)dX m '. If the tangent vector v is to be Cauchy 
characteristic these forms must vanish. This leads to the condition 3.4. The 
conditions 3.5 now follow from inserting 3.4 into the Cauchy characteristic 
conditions 'v J de aug = O' and using 3.6. 

q.e.d. 

Remark. This result expands the treatment to the ’constrained’ case that 
Cartan gave for the ’unconstrained’ variational problem in "Lecons sur les 
Invariants Integraux". See [2] for the connection with the traditional 
'Lagrange Variational Problem' as expounded in Caratheodory's book and 
for the definition and properties of 'symplecdc foliations' and further 
detail. 

4. The differential equations of constrained mechanics. 

There is considerable confusion in the Literature beween the 
Lagrange Variational Problem (or 'constrained extrema') and 'constrained' 
(and 'non-holonomic') mechanics'. I will now describe the latter. Suppose 
again given the following data: 


A smooth paracompact manifold X 


(4.1) 
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T(X) is its tangent vector bundle (4.2) 

A set {0, col, a>m} of smooth 1-forms on X. (4.3) 

Definition. Let {x: t --> x(t)} be a curve in X. It is said to be a trajectory 
of the constrained mechanical system associated with the data 
4. 1-4.3 iff. the following conditions are satisfied: 


0 = o>l(dx/dt) = ... = co m (dx/dt) (4.4) 

There is a curve in X x R m of the form 
(t --> (x(t), fii(t), ..., Mm(t)} such that: 

[dx/dt] J de = m(t)coi + _ + ^(t) com (4.5) 

In other words, 4.4-4. 5 define an ODE system whose solutions are curves in 
X x R m . The 'constrained mechanics trajectories' are the projections in X 
under the Cartesian map projection {X x R m — > X} of the solution curves of 
the ODE system 4.4-4.5. 

5. 'Holonomic' constraints. Equivalence of the Constrained 
Extremal and Mechanics equations in the 'holonomic' case. 

Suppose given the following data: 

A smooth paracompact manifold X (5.1) 

A set {0, co 1 , ..., com} of smooth 1-forms on X. (5.2) 

Indices 1 ^ a, b, ... s m (5.3) 

Definition. The constraint forms {co a } are said to be holonomic iff. there 
is a matrix }co a b} of 1 -forms such that: 


dco a = -2-b co a b A ^ 


(5.4) 
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Remark. Locally, condition 5.4 is equivalent to the following, more 
'geometric', condition: 

The Pfaffian System {w a = 0} is Frobenius Integrable (5.5) 

Let us now combine conditions 5.4 and the Constrained Extremal 
equations 3.5. The following equations result: 

dx/ dtj d0 = - 2 a b^a(t)[dx/ dt] J (<o a bA tob) - 2 a [dx a (t)/ dt] u> a (5.6) 

Rewrite this as follows: 

dx/dtj de = - 2abXa(t)[(u) a b(dx/dt) cob) . a)b(dx/dt) co a b (5.7) 

■ 2]a[dXa(t)/ dt] co a 

The second term on the right hand side of 5.7 vanishes as a consequence of 
the Constraint Equations 4.4, resulting in the following: 

[dx/dt] J de - - 2 a iAa(t)[(co a b(dx/dt) cob) . 2 a [dx a (t)/dt] co a (5.7) 

Theorem 5.1. Let 5.4 be satisfied and let the curve {t --> x(t)| be a 
solution of the Constrained Extremal Equations. Then, { t — > x( t) } is also a 
solution of the Constrained Mechanics Equations 4.4-4. 5. 

Proof. That functions {t --> n a (t)] exist satisfying 4.5 is evident from 5.7. 

q.e.d. 


Here is the converse: 

Theorem 5.2. Let 5.4 be satisfied and let the curve {t --> x( t) } be a 
solution of the Constrained Mechanics Equations 4.4-4.5. Then, {t — > x(t)| 
is also a solution of the Constrained Extremal Equations. 
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Proof. We must show that the existence of functions {t --> X a (t)} satisfying 
5.7 is a consequence of the existence of functions ft --> y. a (t)} satisfying 
4.4-4.5. Examining the right hand side of 5.7, we see that the fn a (t)} can be 
obtained by solving an ODE whose coefficients depend on the {X a (t)}. 

q.e.d. 

6. The constrained mechanics equations in a 'Hamiltonian' form. 

So far, we have been working in the context of general manifold 
theory. Let us specialize now to the situation which is close to the 
’Hamiltonian’ formalism in the traditional particle mechanics case. 


Suppose given the following data: 

n is an integer (6.1) 

The following range of indices: (6.2) 

1 s i, j, ... sn 

X = R 2n+ 1 = R n x Rn x R (6.3) 

{qi, Pi, t} are Cartesian coordinates on X. (6.4) 

{(q, p, t) --> H(q, p, t)} is a smooth real-valued (6.5) 

function on X, called the Hamiltonian. 

0 = 2ipidqi - Hdt (6.6) 

dH = 2iHidqi + 2iHidpi + H t dt, (6.7) 


where {Hi, H*, H t } are the partial derivatives' of 
the Hamiltonian function with respect to the 
'canonical' coordinates 6.4. 
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Theorem 6.1. d0 = 2i(dqi- Hklt)A(dpi + Hidt) (6.8) 

Proof. Follows from 6.7 and 6.6, by a direct compilation, which is left to te 
reader. 

Theorem 6.2. Let V be a smooth vector field on X. Then: 

VJ de = Zi(V(qi) - H*V(t))(dpi + Hidt) - 2i(V( Pi ) + HiV(t))(dqi - Hkit) (6.9) 
In particular, if: 

V(t) - 1 (6.10) 

then: 

VJ de = 2i(V(qi) - Hi)(dpi + Hidt) - 2i(V(pi) + Hi)(dqi - Hkit) (6.11) 

Proof. Apply the operation 'VJ ' to both sides of 6.8. 6.11 follows from 
substituting 6.10 into 6.9. 

Theorem 6.3. Keep the hyopotheses of Theorem 6.2 and condition 6.10. 
Suppose further that: 

VJ de = nUjSiajdq 1 (6.12) 

where {afi are smooth functions on X and ft --> n(t)J is a real-valued 
function of 't'. Then, the following relations must be satisfied: 


V(qi) = H 1 

(6.13) 

V(pO + Hi - n(t)3i 

(6.14) 

Zi(V(pi) + Hi)Hi = 0 

(6.15. 



01/10/95 12:47 ©617 258 8682 


LCS/AI MIT EDU 




TTU MATH DEPT 


@013 


45 Constrained Mechanics and Variational Problems 

Proof. 6.13-.15 results from combining 6.11 and 6.12, and comparing 
coefficients of independent coordinate differentials on both sides of the 
resulting differential form relation. 

Theorem 6.4. Let V be the vector field on X defined by 6.10 and 6.12. 
Then, the orbit curves {t --> (q(t), p(t), t)} of V are solutions of the 
following ODE's: 


dqVdt = dH/dpi 

(6.16) 

dpi/dt = - dH/dq 1 + n(t)ai 

(6.17) 

2iai(p(t), q(t), t)[dqi/dt] = 0 

(6.18) 


Proof. Follows from 6.13-6.15. 

Remark. Equations 6.16-6.18 form an ODE system of (2n+l) equations for 
the (2n+l) 'unknowns: {pi(t), qi(t), |i(t)}. They are the Hamiltonian 
version of the Lagrange Equations of Motion for Constrained 
Mechanics. (In this case, there is only one 'constraint, namely 6.18. The 
case of more constraints can be handled similiarly.) 

7. Final remarks about generalizations. 

The material in Section 6 suggests a Generalization of material about 
Symplectic Manifolds, Geometric Quantization, etc. from the traditional case 
abstracted from Particle Mechanics (as in the work of Dirac, van Hove, 
Segal, Kostant, Souriau, etc) to a abstract sitaution paralleling the material 
developed in Section 6. 

I will briefly sketch such generalizations. Instead of the 'R2n+i' 
situation of Section 6, suppose that we are on a manifold X, with the 
following relation: 


de = Q - dHAdt 


(7.1) 
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'IT and 't' are smooth functions on X. 6 and Q are, respectively, a 1- and 2- 
form on X and Q is closed. Suppose that Vh is a vector field on X such that: 

VrJ de = (xco (7.2) 

V H (t) = 1, (7.3) 

where V is a 1-form defining the constraints and V' is a function on X. 

7. 1-7. 3 imply: 

V H J Q - V(H)dt + dH = nco (7.4) 

This relation generalizes the duality relation between 'infinitesimal 
symmetries' and 'conserved functions’ that plays the basic role in the 
'geometric quamtization' theory of unconstrained conservative mechanical 
systems. I plan to study this Geometric Structure further in a later 
Volume. 
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